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ABSTRACT 

We develop a general expression for the y — y absorption coefficient, for y-rays propagating in an arbi- 
trary direction at an arbitrary point in space above an X-ray-emitting accretion disk. The X-ray intensity is 
assumed to vary as a power law in energy and radius between the outer disk radius, Rq^ and the inner radius, 
which is the radius of marginal stability for a Schwarzschild black hole. We use our result for to 
calculate the y — y optical depth, for y-rays created at height z and propagating at angle <D relative to the 
disk axis, and we show that for d) = 0 and z > Rq, oc where a is the X-ray spectral index and E 

is the y-ray energy. As an application, we use our formalism to compute the minimum distance between the 
central black hole and the site of production of the y-rays detected by EGRET during the 1991 June flare of 
3C 279. In order to obtain an upper limit, we assume that all of the X-rays observed contemporaneously by 
Ginga were emitted by the disk. Our results suggest that the observed y-rays may have originated within 
<45GM/c^ from a black hole of mass ^10^ Mq, perhaps in active plasma located above the central funnel 
of the accretion disk. This raises the possibility of establishing a direct connection between the production of 
the observed y-rays and the accretion of material onto the black hole. We also consider the variation of the 
optical depth as a function of the angle of propagation O. Our results indicate that the “focusing” of the 
y-rays along the disk axis due to pair production is strong enough to explain the observed degree of alignment 
in blazar sources. If the y-rays are produced isotropically in y-ray blazars, then these objects should appear as 
bright MeV sources when viewed along off-axis lines of sight. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: active — 

galaxies: individual (3C 279) — gamma rays: theory — radiative transfer 


1. introduction 


Over the past several years, the EGRET instrument on board the Compton Gamma Ray Observatory (CGRO) has detected 
high-energy y-rays from 44 active galactic nuclei (AGNs), and several of these objects have also been detected at lower energies by 
COMPTEL and OSSE (Fichtel et al. 1994; Hartman et al. 1992; Hermsen et al. 1993; McNaron-Brown et al, 1994; von Montigny 
et al. 1995). One of the most dramatic events seen was the intense y-ray flare of 3C 279 observed in 1991 June, which has called into 
question the viability of a number of theoretical models for y-ray emission in AGNs. Most of the energy emitted during the flare 
appeared at energies exceeding ~ 100 MeV, with an implied (isotropic) luminosity of ^ lO'^® ergs s“‘. However, the actual 
luminosity may have been several orders of magnitude lower if the emission was not isotropic, but rather “focused” onto a small 
region of the sky due to some alignment mechanism such as relativistic beaming or geometrical collimation. In particular, the 
hypothesis of relativistic motion is consistent with observations of apparent superluminal motion, flat radio spectra, rapid optical 
variability and high polarization in blazars. 

In addition to beaming and geometrical collimation, y — y pair production in the X-ray field of the accretion disk could in 
principle lead to preferential escape of the y-rays along the symmetry axis of the disk, due to the strong angular dependence of the 
pair-production cross section. The phenomenon of y - y “focusing” is related to the more general issue of y - y transparency, 
which sets a minimum distance between the central black hole and the site of y-ray production, z^^. One usually finds thatz^^ greatly 
exceeds the cooling length. 
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( 1 . 1 ) 


for relativistic electrons of energy E traveling at essentially the speed of light c and suffering inverse-Compton and synchrotron 
losses in the intense magnetic and radiation fields close to the center of the active nucleus (Dermer & Schlickeiser 1993). An efficient 
pair creation and/or reacceleration mechanism is therefore required in order to explain the presence of relativistic electrons 
sufficiently far from the black hole to produce observable y-rays. The nature of this mechanism is unclear, and several possibilities 
have been suggested, including shocks propagating along a jet (Blandford 1993), the decay of charged pions resulting from hadronic 
collisions (Bednarek 1993), and the acceleration of electrons and positrons via wave-particle interactions in the funnel of a thick 
accretion disk (Becker, Kafatos, & Maisack 1994). Constraints on the mechanism are provided by estimates of the distance over 
which it must operate, z^y — Clearly, the closer one can place the y-ray production site to the black hole while maintaining 
y-ray transparency, the easier it will be to power the observed emission using accretion. It is therefore interesting to calculate the 
minimum distance satisfying the y-ray transparency constraint in the context of a specific X-ray emission scenario, and that is the 
focus of this paper. 
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Gould & Schreder (1967) calculated the absorption coefficient (a^y x cm“ describing the process y y' e'^ -h e~ for high- 
energy photons traversing a soft photon gas, and their formalism has been used subsequently by Blandford (1993), Dermer & 
Schlickeiser (1994), and Bednarek (1993) to estimate the absorption probability for y-rays traversing the X-ray field of an AGN, The 
detailed results depend strongly on the energy and angular distribution of the soft radiation and are therefore sensitive to the 
presence of scattered radiation, which is more isotropic than the primary radiation emitted by the disk. Blandford (1993) and 
Dermer & Schlickeiser (1994) assumed that the soft photon distribution is isotropic, which is reasonable if the X-ray distribution is 
dominated by the scattered component. However, the relative importance of the scattered photons compared to the primary 
(accretion disk) photons is determined by the electron scattering optical depth above the disk, which is not known with any 
certainty. Bednarek (1993) adopted the opposite viewpoint in his calculation of the optical depth to pair production (t^^), neglecting 
electron scattering and instead treating the full anisotropy of the soft radiation field produced by the accretion disk. He assumed a 
standard disk radiating as a blackbody with local temperature T{R) x R The resulting X-ray intensity (integrated over the face 
of the disk) has the characteristic power-law shape x (e.g., Pringle 1981). This is inconsistent with most AGN (and blazar) 
X-ray spectra (with x typically), which are usually better fitted using hot, two-temperature accretion disk models (e.g., 
Shapiro, Lightman, & Eardley 1976, hereafter SLE) in the energy range relevant for pair production (€ > 261 eV for y-ray energy 
E - lGeV;seeeq. [3.4]). 

In y-ray loud blazars such as 3C 279, variability arguments suggest that a significant fraction of the observed X-rays are probably 
produced in a relativistic beam moving toward the observer (Makino et al, 1989). Since these X-rays are co-aligned with the 
observed high-energy y-rays, they do not contribute substantially to the y — y opacity, and therefore they do not strongly constrain 
the minimum distance between the black hole and the site of y-ray production. However, as we demonstrate below, an estimate of 
this distance can be obtained by considering the attenuation of the y-rays due to collisions with X-rays emitted by an underlying 
accretion disk. In order to tie the distance estimate we obtain to the multifrequency observations of particular objects such as 
3C 279 without creating unnecessary theoretical complexity, in this paper we make the further simplifying assumption that all of the 
contemporaneously observed X-rays were emitted by the accretion disk. By overestimating the X-ray flux emitted by the disk, we 
obtain an “ upper limit ” on the minimum distance between the black hole and the point of creation of the observed y-rays, and we 
expect that a proper accounting for the component of the X-ray flux generated in the relativistic beam would further reduce this 
distance. In any event, the upper limit we obtain is small enough to allow for the possibility of establishing a direct connection 
between the production of the observed high-energy y-rays and the accretion of matter onto the black hole. 


2. X-RAY INTENSITY OF THE DISK 

Electrons in the hot, inner region of a two-temperature accretion disk produce power-law inverse-Compton X-ray emission by 
upscattering UV radiation produced in the cool (single-temperature) surrounding region. The temperature of the electrons in the 
hot region is maintained at a nearly uniform value in the range ^ 10^-10^® K by the “Compton thermostat” mechanism (SLE; 
Eilek & Kafatos 1983, hereafter EK), and the emitted X-ray spectrum has the characteristic power-law shape 


L X 6 



(2.1) 


in the frame of the source, where y = (4kTJm^c^) Max (r^j., ^ 1 is the Compton y-parameter, is the electron mass, and is 

the electron-scattering half-thickness of the disk. The power law extends up to photon energies e ~ k7^, above which the spectrum 
declines exponentially. Since a and are essentially constant in the hot inner region of the disk, it follows that the radial variation of 
the X-ray intensity follows that of the total flux, given by 


F{R) = 


3 GMM r _ /6GMY^^' 
871 R^ I ~\ ) 


( 2 . 2 ) 


for a quasi-Keplerian accretion disk around a Schwarzschild black hole of mass M, where Af is the accretion rate (Shakura & 
Sunyaev 1973). Note that the total flux has a nearly power-law dependence on radius F(R) x \ with a sharp cutoff at the radius 
of marginal stability for circular orbits, R^^s = 6Rg, where Rg = GMjc^. 

Based upon equations (2.1) and (2.2), we conclude that the X-ray intensity in the frame of the disk can be represented approx- 
imately by 


/,(€, R) = (oc ergs s ‘ cm ^ sr ' ergs ') , (2.3) 

where a is given by equation (2.1), Rq is the outer radius of the two-temperature region, and R^^ < R < Rq. According to equation 
(2.2), o> = 3 for a two-temperature disk, but we shall leave co as a free parameter so that we can make comparisons with uniformly 
bright disks, with w = 0. Since it is unclear to what degree electron scattering isotropizes the X-rays above the disk, we follow 
Bednarek’s approach, and neglect the effect of electron scattering on the angular distribution of the X-rays. Our neglect of a 
scattered X-ray component is in some sense self-consistent, since it turns out that most of the y — y attenuation occurs very close to 
the disk, where the X-ray distribution is thought to be dominated by unscattered accretion disk photons (Dermer & Schlickeiser 
1994). However, the scattered radiation may attenuate the y-rays over much larger length scales. 

Although the inner region of the disk terminates at the radius of marginal stability, the outer radius Rq of the two-temperature 
region is more difficult to quantify because of its dependence on the value (and possible variation) of the viscosity parameter a, but 
one typically finds 30K^ < Rq < \00Rg (SLE; EK). Most of the photons emitted in the “cool” region beyond radius Rq soft 
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to produce pairs via collisions with y-rays of energy E ^ 1 GeV, and therefore the effect of these photons on the y — y absorption 
coefficient is negligible. 

In order to link our calculation of the y — y absorption coefficient with the X-ray and y-ray data, we must express the disk-frame 
intensity parameter Iq in equation (2.3) in terms of the parameter Fq characterizing the observed (cosmologically Doppler-shifted) 
X-ray flux. 




(oc s * cm ^) , 


which is related to the observed intensity 


— 2) b" ’ 


(2.4) 


(2.5) 


by 

f; = I /; cos 0'dii' , (2.6) 

where primes denote quantities measured in the frame of the Earth. 

Since we are mainly interested in blazar sources, where we are presumably viewing the accretion disks face-on, we can express Fq 
in terms of Iq and the luminosity distance to the source D by substituting for the intensity in equation (2.6) using equation (2.5) and 
integrating over the face of the disk. Making the substitutions dO! = 2n sin O'dO' and tan 6' = R/D yields 


or, in the limit D P Rq, 



(2.7) 


( 2 . 8 ) 


for oj # 2, where Iq is related to I q via the Lorentz transformation of the intensity, which gives for a source located at redshift (e.g., 
Begelman, Blandford, & Rees 1984) 

/o = (l +ir)-^-“/o- (2.9) 

In § 3 we obtain a general formula for the y — y absorption coefficient for y-rays propagating above an accretion disk radiating with 
the X-ray intensity distribution given by equation (2.3). 


3. PHOTON-PHOTON PAIR PRODUCTION 

The exact cross section for y — y pair production is given by the well-known result (e.g., Jauch & Rohrlich 1955), 


<^yy{P) = '^t(1 - P^) 


(3-r)ln(^|^j-2/?(2-j?^) 


(3.1) 


where P is the velocity of the positron (or electron) in the center-of-momentum (CM) frame in units of c, and <Tj is the Thomson cross 
section. The Lorentz factor of the particles in the CM frame, T, is related to the energy of the soft photon e, the energy of the y-ray £, 
and the interaction angle between the photon trajectories 0 by 


1 eF(l-cos6>) 
1 — 


(3.2) 


where 6, F, and 6 are all measured in the same (arbitrary) frame. The cross section for the process, (7^^,, vanishes as T 1. Hence for 
fixed e and £, the absorption probability decreases as 0-^0, and consequently the optical depth to pair production (t.,^) drops 
sharply as the photons become co-aligned. We therefore anticipate that the value of t,,y will depend sensitively on the angular 
distribution of the soft radiation relative to the direction of propagation of the y-ray. 


3.1. y — y Absorption Coefficient 

Working in terms of the specific intensity f, we can write the differential contribution to the absorption coefficient due to soft 
photons propagating in energy range de and in directional range dQ as 
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We assume that the X-ray intensity of the disk is given by equation (2.3), and that any change in the X-ray spectral index a occurs 
below the minimum soft-photon energy required to produce a pair with F = 1, 




2m^c^ 


£(1 — cos 6) 


(3.4) 


Taking E = 1 GeV, we find that = 0.522 keV for a right-angle collision (cos 6 = 0), and = 0.261 keV for a head-on collision 
(cos 6 = “1), which is the minimum possible X-ray energy resulting in pair production. The X-ray spectrum given by equation (2.3) 
would conflict with the measured EGRET spectrum if extended to arbitrarily high energies, and consequently a break must exist in 
the energy range 0.1-1 GeV. However, for our purposes, the X-ray power law can formally be extended to infinite energy without 
producing numerical divergence because oc € ” Mn 6 as 6 -► oo for fixed E and 0. 

Using equation (2.3) for the intensity, integration of equation (3.3) over e and Q yields for the absorption coefficient 


a 


yy 


Iq 



- cos 6)d€dQ , 


(3.5) 


where 0 is the angle between dQ and the direction of propagation of the y-ray, and the domain for the n integration corresponds to 
the range of solid angles subtended by the disk. Equation (3.5) can be broken into two separate integrals by transforming from e to ji 
using equation (3.2). The result is 


lo 

= 7 


where the integration over P is expressed by the function 


a r / R\~^ 

I (Tj ^(a) J (1 - cos ’ 


2p{\ ~ p^r-^<T^,(p)dp . 


(3.6) 


(3.7) 


plotted in Figure 1. 

Equation (3.6) can be used to calculate the absorption coefficient for a y-ray propagating in an arbitrary direction at an arbitrary 
point in space. However, it is sufficient for our purposes to consider y-rays propagating in a plane that includes the symmetry axis of 
the disk. Referring to the geometry indicated schematically in Figure 2, we assume that the y-ray propagates in the (x, z) plane, 
where the z-axis is the symmetry axis of the disk, and we set the azimuthal angle <p = 0 along the x-axis. In general, we wish to relate 
dQ to the differential area in the disk dA = RdR d(p using 


cos tj RdRd(p 
dCl = 


(3.8) 


where rj is the angle between the X-ray trajectory and the normal to the disk surface, and / is the distance between the photon- 
photon interaction point and the X-ray emission point in the disk. Representing the location of the y-ray using polar angle S and 
distance from the origin p, and assuming that the y-ray trajectory makes an angle <I> with the z-axis (see Fig. 2), our expression for the 



0 12 3 4 


a 

Fig. 1. — The cross section integral function 4^(a)(eq. [3.7]). For the 3C 279 data, o: = 0.68 (see Makinoet al. 1993), and 'P = 0.245. 
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X-ray trajectory f 



Fig. 2 . — Schematic representation of the 7-ray propagation geometry. The y-ray is located at distance p from the origin and polar angle d and propagates in a 
plane that includes the 2-axis. Attenuation of the y-rays occurs as a result of collisions with X-rays generated at all points on the disk. The interaction angle between 
the photons is 0, and the angle between the y-ray trajectory and the z-axis is <I>. 


absorption coefficient becomes 


where 


cc JE, p, 3, 0)) = 


In 


2m ^ c‘ 






(1 


cos rjdipRdR 

Uj — 


/2 = ^2 _ 2pfi sin ^ cos (p , 


p cos (3 — <t>) ~ R cos (p sin 0 
cos 6 — 


p cos 3 
cos ri = — 


(3.9) 


(3.10) 


Note that oc which verifies that softer X-ray spectra yield larger optical depths. Equations (3.9) and (3.10) comprise our 
fundamental result for the y — y absorption coefficient for a y-ray of energy E propagating in the power-law X-ray field of an 
accretion disk. 


3.2. Absorption Coefficient for Propagation along the Disk Axis 

One of the most interesting cases to consider is that of the attenuation of y-rays propagating along the symmetry axis of the disk, 
since this scenario is of direct relevance for the interpretation of the y-rays detected by EGRET from blazar sources. In this case 
equation (3.9) is simplified substantially due to the cylindrical symmetry, and the integration over <p is trivial. We can obtain the 
absorption coefficient for a y-ray of energy E propagating outward along the symmetry axis at height z by setting p = z and 
^ = (I) = 0 in equation (3.9) and transforming from Rio 6 using tan 6 = R/z, yielding 

oCy/E, z) = J ^ tan“" 9 sin 6d0 , (3.1 1) 

where 0^ax(^) [^min(^)] the maximum (minimum) value of the interaction angle 0 at height z. Neglecting gravitational curvature of 
the photon trajectories, we have 


By writing equation (3.1 1), we have implicitly assumed that z > h, where h is the half-thickness of the disk. This condition is trivially 
satisfied for a thin disk, but the disk may also be considered to be thick, so long as we confine our attention to points sufficiently far 
above the surface. 
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4. CALCULATION OF THE OPTICAL DEPTH 

We can use equations (3.9) and (3.10) to calculate the y - y optical depth for a y-ray of energy E created at radius R and height 
z and propagating at angle <I> with respect to the z-axis by writing 


'*0C’ 

Tyy(£, R, Z, = (Xyy(£, p, S, <l>)dA , 

Jo 


where A is the distance traversed by the y-ray since its creation, and p and <5 are given by 

, , , ^ , R A sin <t> 

^ ,2 ^ 2;.(^ sin O + z cos d>) , sin ^ ■ 


cos 3 = 


+ A cos CD 


P P 

In keeping with equation (3.9), we have assumed that the plane of propagation of the y-ray includes the z-axis. 


(4.1) 


(4.2) 


4.1. Propagation along the Disk Axis 

We can obtain some further simplification for y-rays propagating along the z-axis, since then 3 = = R = 0, and p = z The 

optical depth to pair production for a y-ray of energy E created at height z above the center of the disk and propagating outward 
along the z-axis is then given by 

Ty.,(£, -) = 

where is given by equation (3.1 1). Unfortunately, one must resort to multidimensional numerical integration to compute even 
in this highly symmetrical case. However, it is possible to obtain a closed-form expression for Tyy if z > Rq. Approximating the 
angular integral in equation (3.1 1) yields for the absorption coefficient 


ayy(£, z) = A 



(4.4) 


where 


A = 


nl o(Tt 4"(a) 


^ 2flr + 4 - tt> 


t) ] 


(2a -I- 4 — cu)r 

Substituting equation (4.4) into equation (4.3) then yields the approximate far-field solution for the optical depth 


Tyy(£, Z) ^ 


2a + 3 


(riij.)' 


(4.5) 


(4.6) 


valid if z > - In § 5 we use this result to estimate the height of the y-ray " photosphere,” defined as the value of z such that Tyy = 1. 

We also calculate the y — y optical depth for the y-ray blazar 3C 279 by focusing on three illustrative scenarios for the propagation 
of the y-ray: (i) propagation along the disk axis; (ii) propagation along a line intersecting the disk axis; and (iii) propagation parallel 
to the disk axis. 


5. APPLICATION TO 3C 279 

The observations of strong fluxes of X-rays and high-energy y-rays from the blazar 3C 279 in 1991 June allows us to place 
constraints on the location of the y-ray emission site based on y ~ y transparency requirements. Under the assumption that the 
contemporaneously observed X-rays were produced by a two-temperature accretion disk, we can use the expressions developed in 
§§ 3 and 4 to calculate the optical depth traversed by the observed high-energy y-rays. Since variability arguments (Makino et al. 
1989) suggest that a significant fraction of the X-rays were probably produced in a relativistic jet, the results we will obtain for the 
optical depth are upper limits. 

The X-ray flux reported by Makino et al. (1993; Fig. 2) during the outburst had spectral index a = 0.68 and flux parameter 
Fq = 5.76 X 10“^ s'^ cm^^ (note that the flux given by their eq. [1] appears too large by a factor of 10). For this value of a, we 
obtain y = 1.6 (see eq. [2.1]) and 'F(a) = 0.245 (see Fig. 1), and for the redshift of 3C 279 = 0.538), we obtain the luminosity 

distance (Lang 1980) D = Hq + 0.5(1 - go)^^] = 5-65 x 10^“^ cm assuming = 0.5 and Hq = 100 km s “ ‘ Mpc L We can 
relate Fq to the disk-frame intensity parameter Iq using 

according to equations (2.8) and (2.9). 


5.1 . Distance to the Gamma-Ray Emission Site 

As an application, we consider first the simplest geometrical case, that of a y-ray propagating away from the black hole along the 
axis of a two-temperature disk emitting with the power-law X-ray intensity given by equation (2.3). This scenario is useful for 
determining the minimum height above the disk at which the observed y-radiation could have been produced in 3C 279 in 1991 
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June. We can calculate the optical depth to pair production for a 7-ray of energy E created at height 2 above the center of the disk 
and propagating outward along the z-axis using equation (4.3). In order to calculate we must assume values for the inner 
and outer (Rq) radii of the two-temperature, X-ray-emitting region. In Figure 3 we set = 6R^ (accretion onto a Schwarzschild 
black hole) and Rq = 30 Rg (a typical value for the outer radius; see SLE and EK), and plot log as a function of log (M/Mq) and 
log (z/Rg) for a 7-ray of energy E = 1 GeV using the X-ray data taken during the 1991 June EGRET flare of 3C 279. The optical 
depth exceeds unity in the cross-hatched area, and = 1 along the heavy solid line denoting the y-ray photosphere, with height 
2 = Zyy, For reference, 2 = Rq along the upper dashed horizontal line, and 2 = R^, along the lower dashed line. The slope of the 
contours depends on how close to the disk the y-rays are produced relative to R^^ and Rq . For 2 < R^^, the contours become 
increasingly vertical because the absorption coefficient loses its dependence on 2 close to the disk. For 2 > Rq, the contours make a 
transition to the far-field power-law behavior given by equation (4.6), which reduces to Ty,, = 2.4 x 10^^(M/Af q)“ ^(z/R^)""^ for 
the parameters used in Figure 3. 

In the case of 3C 279, we find that when Af > 4 x 10‘^ Mq, the intensity is so low that < 1 regardless of the value of 2. On the 
other hand, we see that observable 7-rays can be produced quite close to the black hole for reasonable values of M if 2 > Rq. For 
example, if M = 10^ Mq, then <1 for z > z^^ = 45R^ = 2 x 10 pc. This suggests that y-rays produced at the top of the 
evacuated funnel at the center of a hot accretion disk (perhaps as a consequence of MHD waves generated in the surrounding 
hot-disk plasma; see Becker & Kafatos 1994) can escape without y — y attenuation by the disk X-rays (although they may suffer 
attenuation by scattered X-rays on larger spatial scales; see Dermer & Schlickeiser 1994). 

The precise value of the outer radius Rq of the two-temperature region is uncertain, but it is likely to fall in the range 
30 < Rq/R^ ^ 100 (SLE; EK). In Figure 4 we therefore recompute the contours using the same parameters as in Figure 3, except 
Rq = lOOR^. Here again 2 = Rq and 2 = R ^^^ along the upper and lower dashed lines, respectively. The differences between Figures 3 
and 4 are not dramatic, although the height of the y-ray photosphere has increased from Zy^ = 45Rg to Zyy = lORg for Af = 10^ A/q. 
This is due to the presence of X-ray photons emitted at large radii, for which the interaction angle 0(and therefore the probability of 
pair production) is enhanced. In this case, equation (4.6) for the approximate far-field solution reduces to ^ 3.6 
X 10^^(AJ/AfQ)“^(z/R^)''^ in good agreement with Figure4forz > Rq. ^ 

To investigate the effect of the angular dependence of the y - 7 absorption coefficient on the values obtained for in Figure 5 
we plot the height of the y-ray photosphere, z^/R^, as a function of the black hole mass M/Mq for R^, = 6R^ and a; = 0 (uniformly 
bright disk; dashed line) or to = 3 (two-temperature disk; solid line). We compare results obtained using our two standard values for 
the outer radius, Rq = 30R^ (open squares) and Rq = lOOR^ (open triangles). When cu = 3, increasing Rq from 30R^ to lOOR^ causes 
only a modest increase in the height of the photosphere as mentioned above. However, going from cu = 3 to cu = 0 while holding Rq 
constant results in a substantial increase in Zyy due to the increased amount of X-radiation impinging on the axially propagating 
y-rays from large radii, where the interaction angle 0 has its maximum value. This effect is particularly strong for Rq = lOOR . 

It is also interesting to consider the effect of varying the inner radius of the disk, R^^, since this will also impact on the^^angular 
distribution of the X-rays. In Figure 6 we plot the photosphere height zJRg as a function of M/M q for R^JRg = 2, Rq/R^ = 
(30, 100), and a> = (0, 3). These results may represent approximately situations involving disks around Kerr black holes, for which 
the radius of marginal stability can approach 1.2R^ (Shapiro & Teukolsky 1983). In this case equation (2.2) for the flux F(R) is no 




Fig. 3. — Contour plot of the log of the y - y optical depth (eq. [4.3]) for y-rays with energy E = I GeV created at height c above the center of a 

two-temperature disk (<y = 3) with outer radius = 30/?^ and inner radius = 6R^. The upper and lower dashed horizontal lines denote z = /?o and z = R^^, 
respectively, and the y-rays propagate outward along the z-axis. The curves were computed using the X-ray data for 3C 279 taken during the 1991 June EGRET flare 
by Makino et a). (1993). 

Fig. 4. — Same as Fig. 3, except Rq = lOOR^ 
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log(M/3/g) \og{M/M(^ 

Fig. 5. — Height of the y-ray photosphere, plotted as a function of the black hole mass MjM^ for = 6R^ and to = 0 {dashed line), a) = 3 (solid line), 

Rq = 30Rg {open squares), Rq = \00R^{open triangles). The curves were computed using the 1991 June 3C 279 X-ray data. 

Fig. 6. — Same as Fig. 5, except = 2R^ 


longer valid, although the flux still maintains the approximate power-law radial dependence F(R) oc R ^ (Eilek 1980). Obviously 
our model cannot be rigorously defended in this context since we have completely neglected relativistic effects, but, bearing this 
caveat in mind, it is interesting to note the modest decrease in the height of the photosphere going from Figure 5 to Figure 6 when 
^ ^ 3 ^ to the decrease in the flux-averaged value of the interaction angle 6. By contrast, the results depicted in Figures 5 and 6 
agree closely when cu = 0 (uniformly bright disk), because then most of the emission comes from radii R > R^^, and the value of 

is therefore irrelevant. . i- 

In Table 1 we give values for the height of the y-ray photosphere Zyy/Rg, along with values for the associated light-crossing 
timescale (which may approximate the variability timescale), f,, - zjc, for £ = 1 GeV, w = 3, and various values of the black hole 
mass M/M Qy the outer disk radius Rq/R^, and the inner disk radius R^^/R^. Table 1 also includes values for the disk-frame intensity 
parameter Iq. Using our standard parameters M = 10^ M©, Rq = 30Rg, R^^ = 6Rg, and cu = 3, we find that ^ 2 days, in 
agreement with the y-ray variability timescale reported by Kniffen et al. (1993). 

5.2. Gamma-Ray Collimation 

The dependence of the y - y absorption probability on the angular distribution of the X-rays suggests that pair production may 
help to collimate the high-energy radiation escaping from blazars. We expect that for given values of the inner and outer disk radii 
R„, and Rq, the optical depth to pair production for y-rays originating at a fixed location above the center of the disk will increase 
with increasing propagation angle O, mainly due to an increase in the flux-averaged value of the interaction angle 6, and to a lesser 
extent because the y-rays remain close to the disk for a longer period of time. We shall continue to work with X-ray data obtained 
during the 1991 June flare of 3C 279 in all of the examples. 


TABLE 1 

Model Results for 3C 279 


M 

(Mo) 

^0 

(R.) 

(R,) 

(R,) 

10* 

30 

6 

81.2 


30 

2 

59.9 


100 

6 

141.2 


100 

2 

101.8 

10’ 

30 

6 

45.2 


30 

2 

32.0 


100 

6 

69.9 


100 

2 

45.9 

10“> 

30 

6 

22.5 


30 

2 

14.6 


100 

6 

27.6 


100 

2 

16.3 


t,, A(D /o 

(d) (deg) (s ' cm ^ sr ‘) 


0.46 

5.7 

4.9 X 10^5 

0.34 

6.0 

1.4 X 10^’ 

0.81 

7.8 

4.2 X 10“ 

0.58 

7.7 

1.3 X 10“ 

2.57 

9.5 

4.9 X 10“ 

1.82 

9.7 

1.4 X 10“ 

3.98 

12.1 

4.2 X 10“ 

2.62 

11.6 

1.3 X 10“ 

12.85 

15.3 

4.9 X 10“ 

8.31 

14.6 

1.4 X 10“ 

15.74 

17.2 

4.2 X 10“ 

9.31 

15.4 

1.3 X 10“ 
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We can use equation (4.1) to calculate the optical depth for y-rays of energy E created at a height z above the center of the disk 
and propagating at an angle <t> with respect to the z-axis. In order to explore the angular dependence of y — y attenuation, in Figures 
7 and 8 we plot the transmission coefficient exp ( — calculated using the 3C 279 X-ray data as a function of O for = 6R^, 
E = 1 GeV, M — 10^ Mq, and various values of Rq and z. The transmission coefficient gives the fraction of the y-rays that would 
escape to infinity in a given direction out of those emitted in that direction. We set Rq = ^^R^ in Figure 7, Rq = lOOR^ in Figure 8, 
and, to illustrate the importance of angular effects, we set ca = 3 in Figures la and 8a, and cu = 0 in Figures lb and 8h. Strong 
“focusing” of the escaping y-rays along the z-axis is apparent in all cases, and, for w = 3, the half-width of the angular distribution, 
A<I), falls in the range 10° < AO < 20°. Similar results are obtained for AO when (u = 0 (Figs, lb and 8h), but the overall level of 
attenuation is much greater in these cases due to the enhancement of the X-ray emission produced at large radii in the disk. In Table 
1 we give values for AO evaluated at the photosphere, z = z^^, for various model parameters. Our results for AO are comparable to 
the observationally determined values (e.g., Padovani & Urry 1992), and we therefore conclude that for acceptable disk parameters, 
y — y attenuation alone is enough to produce the degree of y-ray collimation implied by observations which associate extragalactic 
y-ray emission with blazars. 

5.3. Shape of the Gamma-Ray Photosphere 

The high-energy y-rays observed from many blazars are emitted nearly parallel to the axis of the jet, as is evidenced particularly 
by the common occurrence of apparent superluminal motion in the sources. It is therefore interesting to calculate, within the context 
of our model, the y — y optical depth for a y-ray created at radius R and height z, and subsequently propagating parallel to the z-axis 
of the disk, which can be obtained by setting <I> = 0 in equation (4.1). It is particularly interesting to determine the shape of the y-ray 
photosphere by varying the value of the y-ray creation radius R. 

In Figure 9, we plot the photosphere height z^fR^ calculated using the 3C 279 X-ray data as a function ofR/R^ for various values 
of M/Mq, assuming our standard parameters Rq = 30R^, R„, = 6R^, and = 3. We find that the height of the photosphere 
increases dramatically as M decreases, and the shape of the photosphere becomes quite concave. Conversely, for large values of M, 
the photosphere becomes nearly flat. We interpret this behavior as follows. For large values of M, the disk-frame intensity 
parameter /q is small, and therefore the photosphere is located relatively close to the disk, with z^^/Rq < 1. In this case, most of the 
y — y absorption takes place so close to the disk that the flux-averaged value of the interaction angle 6 is nearly 90° at all points in 
the photosphere, and therefore the height of the photosphere varies little across the disk. On the other hand, as M decreases, Iq 
becomes large, and consequently the photosphere lies farther from the disk, with z^fR^ > 1. In this case, the flux-averaged value of 0 
increases significantly going from the center of the photosphere to the edge, and therefore the height of the photosphere must 
increase rapidly in order to keep = 1, as seen in Figure 9. 

6. DISCUSSION 

Although the physics of particle creation and acceleration far from the turbulent plasma accreting onto a black hole is poorly 
understood, there are a number of relatively well-understood plasma and accretion instabilities that can give rise to efficient particle 



$ $ 

Fig. 7. — Transmission coefficient exp plotted as a function of the propagation angle using the 199! June 3C 279 X-ray data for M = 10^ Mq, Rq = 

30Kg, = 6/?^ and (a) <0 = 3, (h) a) = 0. The y-rays are created at the height z/R^ indicated for each curve. 
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acceleration within the disk itself, or within the evacuated funnel region at the center of the disk. For example, resonant turbulent 
acceleration via interactions between particles and (linear or nonlinear) Alfven waves can in principle lead to the production of 
electrons with very high Lorentz factors, and to the associated emission of high-energy inverse-Compton y-rays (e.g., Becker & 
Kafatos 1994; Becker et al. 1994). However, in order for any of these mechanisms to create observable y-rays, the optical depth to 
pair production on the background UV and X-ray photons must be small enough for the y-rays to escape from the source. 

In this paper, we have attempted to bridge the gap between theory and observation by developing useful constraints based on the 
simple hypothesis that the power-law X-ray spectrum originates in the inner region of a two-temperature accretion disk. Such disks 
commonly provide good fits to the hard X-ray spectra of AGNs. In § 3 we developed a general expression for the y — y absorption 
coefficient for y-rays propagating above a two-temperature disk in a plane that includes the symmetry axis, and in § 4 we derived 
an approximate far-field solution for the y — y optical depth for y-rays created above the center of the disk and propagating 
outward along the disk axis, taking full account of the angular dependence of the X-ray distribution. In § 5 we applied our results to 
3C 279 using the X-ray data reported by Makino et al. (1993) during the 1991 June y-ray flare, and demonstrated that the observed 
y-rays could have been produced within <45GM/c^ of a black hole of mass M > 10^ Mq, This is close enough to the central source 
to establish a direct connection between the production of the y-rays and the accretion of material onto the black hole. In particular, 
our results suggest that the observed y-rays may have been produced near the top of the evacuated funnel at the center of a hot, 
two-temperature accretion disk. In this scenario, the relativistic electrons that Compton-upscatter UV photons to produce the 
y-rays are themselves produced as a consequence of photon-photon collisions in the funnel. The absence of protons in the funnel due 
to the centrifugal barrier allows the presence of Alfven waves with wavelengths short enough to effectively accelerate the electrons 
(Becker et al. 1994). 

In order to explore the observational implications of the angular dependence of y — y pair production, in § 5 we also computed z^y 
for y-rays born along the disk axis but propagating at an angle O with respect to it. The results for the half-angle (Ad) < 15°) 
indicated in Figures 7 and 8 and in Table 1 suggest that pair production in the X-ray field of the accretion disk may be strong 
enough to explain the observed degree of y-ray alignment in blazars. Furthermore, the variability (light-crossing) timescale tyy = 
Zyy/c we derive for black hole mass M — 10^ Mo is in good agreement with the 2-day y-ray variability of 3C 279 reported by Kniffen 
et al.(1993). 

The results presented in §§ 4 and 5 suggest that y - y pair production can help to “focus’’ an intrinsically isotropic source of 
y-rays so that the high-energy radiation escapes preferentially along the disk axis. Whether or not the y-rays actually are produced 
isotropically is a separate question. If they are, then several additional considerations come into play, regarding (i) the self-opacity of 
the y-rays; (ii) the fate of the y-ray energy emitted off-axis; and (iii) the required total luminosity. 


6. 1 . y — y Self‘Opacity 

Whether or not the high-energy y-rays are produced by an intrinsically isotropic source, their self-interaction optical depth 
cannot exceed unity if they are to be observed. The most severe self-interaction constraint is obtained in the case of isotropic 
emission. Assuming a monochromatic, isotropic distribution of y-rays with energy Eq > and intensity /, = BS(E — Eq), 
integration of equation (3.3) yields for the self-interaction absorption coefficient 


^self 


where cr.,y(/?) is given by equation (3.1) and 



( 6 . 1 ) 


(6.2) 


If the y-ray source is an isotropically emitting sphere of radius then the intensity parameter B is related to the y-ray luminosity 
Ly via 


B = 


(2nRJ^ 


and we obtain for the self-interaction y — y optical depth 


(6.3) 


^self — 


tiR^cE 


E V* C"- ^ 2lidfi ~ 
.'WfCV Jo <Tt (1 - ' 


(6.4) 


An isotropic y-ray source can be ruled out if z^^^f exceeds unity for reasonable values of R^. The first factor on the right-hand side of 
equation (6.4) is essentially the standard y-ray compactness, and the quantity in brackets becomes 10“^ for £ = 1 GeV. Taking 
Ly = 10'** ergs s“ ^ for 3C 279 and setting < 1 yields the easily satisfied transparency requirement R^ > 10*® cm, and therefore 
self-interaction of the GeV radiation in 3C 279 is not likely to be important. A similar calculation by McNaron-Brown et al. (1995) 
indicates that the self-interaction of radiation in the energy range 0.05-10 MeV (observed by OSSE) is also negligible. 


6.2. Fate of the Off-Axis Radiation 

If the y-rays are produced isotropically, then only those emitted close to the disk axis will escape from the X-ray field of the 
accretion disk. The radiation emitted farther from the axis will be degraded in a pair-Compton cascade (Protheroe, Mastichiadis, & 
Dermer 1992), with most of the luminosity residing ultimately in a broad MeV feature. In this view, objects such as 3C 279 (while in 
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y-ray flare mode) will appear as bright sources of MeV emission when viewed along lines of sight >20“ from the disk axis. Centaurus 
A may be an example of such an object (Skibo, Dermer, & Kinzer 1995). 

The luminosity required to power isotropic y-ray emission, L, = 10** ergs s ' is orders of magnitude higher than that required if 
the emission is beamed. This is perhaps the simplest argument against isotropic emission, but it is especially compelling when 
coupled with the strong circumstantial evidence in favor of relativistic bulk motion. However, there is apparently no strong a priori 
reason to rule out isotropic emission, since attenuation due to collisions with X-rays emitted by the disk is sufficient to produce the 
observed degree of alignment. In any event, our results regarding y - y collimation do not rely on the existence of an isotropic y-ray 
source. The mechanism we consider here can work in concert with other sources of collimation, such as relativistic bulk motion or 
geometrical focusing. 

6.3. Attenuation by Scattered X-Rays 

The backscattering of X-rays produced in the accretion disk by electrons located above the disk tends to isotropize the X-ray 
distribution, and this can have a strong effect on the results obtained for the y — y absorption coefficient. We believe that our neglect 
of a scattered component to the X-ray distribution is reasonable due to the proximity of the y-ray production site to the black hole 
in our model. It has been shown (Dermer & Schlickeiser 1994) that accretion disk photons will dominate over scattered X-rays 
within a distance of ~0.01-0.1 pc from the central source, which is much larger than the typical height of the y-ray photosphere, 
obtained here (see Table 1). However, the problem of y - y attenuation by scattered X-rays on much larger spatial scales is a serious 
issue that must be addressed using detailed calculations that determine self-consistently the angle and energy distributions of both 
the scattered radiation and the electron-positron pairs at large distances from the black hole. Such a calculation is beyond the scope 
of this paper, but it suggests a promising direction for future work. 

Our results indicate that if the emission height z is held fixed, then the optical depth increases with increasing y-ray energy E as 
T X £*, and therefore the height of the y-ray photosphere z,, increases with increasing E. If we arbitrarily set the emission height 
for all of the y-rays equal to the value of z^j. obtained for £ = 1 Ge V, then the observed spectrum will be attenuated by a factor ~ 7 
going from E = 1 GeV to £ = 5 GeV. It is unclear whether such a steepening actually exists in the spectrum of 3C 279 (Hartman et 
al. 1992). Finally, for energies £ > 5 GeV, the interaction with UV radiation produced in the cool outer region of the disk may 
become important, further raising the height of the photosphere. 

To summarize our main points, application of our general expressions describing y - y attenuation above X-ray-emitting 
accretion disks leads us to conclude that in the case of 3C 279, (i) y-ray transparency (to disk-generated X-rays) is achieved if the 
y-rays are produced at height z>Ro ', (ii) y - y attenuation is strong enough to explain the observed association of high-energy 
y-ray emission from AGNs exclusively with blazars; and (iii) the observed y-rays may be produced in active plasma located above 
the funnel at the center of a two-temperature accretion disk. 

We gratefully acknowledge the remarks provided by the anonymous referee, whose insightful comments contributed significantly 
to the paper. This work was partially supported by a NASA Cycle 3 CGRO Guest Investigator grant. 
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